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Abstract Androsace L. consists of more than 100 species distributed mainly in northern temperate 
areas. This genus[] typically characterized by having small flowers with a constricted corolla тоа 
was subdivided into six sections[] sect. Затиећа Schlechtd.[] sect. Mirabiles[] Hand.-Mazz. [| 
Yang & Huang[] sect. Androsace[] sect. Chamaejasme Косћ. [] sect. Aizoidium Hand. -Mazz.[] and 
sect. Orthocaulon Hand.-Mazz. In this genus[] the species with cushion-like growth belong to sect. 
Chamaejasme occurring in alpine habitats. There are other three small genera in the tribe An- 
drosaceae. Pomatosace Maxim.[] with only one species[] P. filicula Maxim.[] endemic to Qinghai- 
Tibetan Plateau[] is traditionally placed in the subtribe Soldanellinae because of its circumscissile 
capsules[] but other morphological characters and the recently reported evidence from cytology indi- 
cate that its affinity is closer with Androsace . Douglasia Lindley has eight homogeneous species in 
Northern American mountains. Vitaliana Sesler[] with only one species in the European Alps[] is 
usually considered a congener of Douglasia . These four general] morphologically similar in having 
small flowers with a constricted corolla throat апі“ Androsace pollen"[] referred here as‘ Androsace 
group’ [] are disjunctly distributed from the Qinghai-Tibetan Plateau[] to Central- and Northern-Asia[] 
Europe and North America. It is an ideal group to elucidate the origin and dispersal patterns of the 
Northern Hemisphere flora[] which has a distribution center in the eastern Qinghai-Tibetan Plateau 
and the Hengduan Mountains. In addition[] the habits of this group range from annuals to rosette 
perennials and cushion-like perennials. The cushion-like species are disjunctly distributed in the 
central Qinghai-Tibetan Plateau and European Alps mountains. Further studies are needed to deter- 
mine whether they are closely related to each other or it is only the result of convergent evolution 





2004-03-100 0 0 2004-08-02 [] 0 000 
gggggpuggggggggggggugaga К$СХ-5%-106ПШПППППППППППППШППППП 
O [L] 30000012[T] 00 Supported by a Knowledge Innovation Project of the Chinese Academy of Science$] Grant. No KSCX-SW- 
106[ T] Foundation for the Author of National Excellent Doctoral Dissertation of Chiná] Grant No. 200327[] and the National Natural 
Science Foundation of Сып] Grant No. 3000001200 

* [] 0 П 00 Author for correspondence. Е-та liujq @ mail. nwipb. ас. en] 





482 uU 00000: 42 





under the habitat pressure selection in separate sites. 

We firstly reported cpDNA ¿rn L-F[] nrDNA ITS sequences respectively or both for 29 popula- 
tions of 14 species in Androsace and Pomatosace mainly from the central Qinghai-Tibetan Plateau. 
A phylogeny was further constructed for these species[] and those distributed in Europe and North 
America in^ Androsace group” and related genera in Primulaceae based on newly reported sequences 
and those downloaded from GenBank. Phylogenetic inferences indicated that the four genera in’ An- 
drosace group” form a well-supported monophyletic clade. Two main clades were discovered in the 

“ Androsace group [] a wholly Androsace clade[] and the other comprising Ротаѓоѕасе] Douglasia 
Vitaliana and nine Androsace species in all analyses. The grouping and position of three species of 
sect. Samuelia and two species of sect. Androsace varied according to the analyses of the different 
datasets[] trn L-F[] ITS or a combination of them. The biogeographical mapping of species distribu- 
tion revealed that the basal species of each clade sited in southeastern China and the eastern Qing- 
hai-Tibetan Plateau[] which obviously supports the origin area hypotheses for" Androsace group” 
based on morphological studies. From the eastern Qinghai-Tibetan Plateau multiple lineages spread 
to and colonized Europe in different ages. The crude timing based on ITS sequence differentiation 
indicated a recent diversification of Androsace group” within Miocene of the Tertiary. The cushion- 
like species evolved independently in Asia and Europe. The origin and diversification of cushion-like 
species in the Qinghai-Tibetan Plateau might be earlier than those in Europe. Their occurrence in 
the Qinghai-Tibetan Plateau area might be correlated with the large-scale uprising of the plateau and 
the mountains building since Miocene. But the appearance of the cushion-like species in Europe 
might result from the fast habitat isolation due to the climate oscillation and the cycles of developing 
and retreating of the ice sheet in the Northern Hemisphere during the Quaternary. The large-scale 
colonization of cushion-like Androsace species in the Qinghai-Tibetan Plateau probably occurred in 
the late Holocene when the global climate re-cooled after the end of glacial cycles. The intraspecific 
genetic differentiation in some species might partly reflect their retreat and re-colonization in the 
Qinghai-Tibetan Plateau area during the climate oscillation. 

Key words Androsace[] Pomatosace[] phylogeny[] biogeography[] Qinghai-Tibetan Plateau[] Heng- 
duan Mountains[] cushion-like plants . 
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01 0000000000 Wendebd] 1961000 00000 00 199000 000 00 0 Џуаћапа [] 
Douglasia 0 D 
Table 1 Source of material$] The infrageneric classification of Androsace is according to Wendelbd] 19610 and Yang 
and Huang] 19900 but excluding Vitaliana and Роигћача 


































































































000 000 00005 000002 свевак000 Ü 
Тахоп Locality Voucher! Distribution ^ GenBank accession number 
+20 
and habit 1$ tmL-F 
0000 Androsace L. 
OOOO Sect. Androsace 
00000 A. filiformis Retz. 0 GenBank[] A/ANN AY275078 AY274955 
00000 A. maxima L. 0 GenBank[] A/ANN AY275099 AY274958 
OOOO A. septentrionalis L. 0 GenBank[] АП A/ANN AY275074 AY274959 
ООО Sect. Aizoidium Hand.-Mazz. 
000000 0 GenBank[] © АММО PRE AY275098 AY274953 
A. integra] Maxim.[] Hand. -Mazz. 
ПОП] 0 Sect. Chamaejasme Koch. 
00000 0000 000 @/ROS AY545758” AY545783 * 
A. alaschanica Maxim. var. zadoensis Nanggén[] Qinghai [] J. Q. Liu[] 585 
Y. C. Yang & R. F. Huang 
00000 00000 000 O/ROS _ AY545733” AY545761 * 
A. bisulca Bur. & Franch. Nagarzé[] Xizang [] J. Q. Liu[] 1125 
0000 000 AY545734” | Ау545762" 
Saneri[] Xizang [] J. Q. Liu[] 1140 
0000 000 AY545735' AY545763 * 
Gyaca[] Xizang [ J. Q. Liu[] 1160 
00000 000000000 000 Q/ROS  AY545756°” AY545779 * 
А. limprichtii Pax & Hoffm. Cultivated in Genevesd] J. Q. Liu[] 708 
Botanical Garden 
00000 0000 000 Q(O/ROS AY545739 AY545764 * 
A. mariae Kanitz Nanggén[] Qinghai [] J. Q. Liu[] 622 
0000 000 AY545737” — AY545765 * 
Qilian[] Qinghai 0 J. Q. Liu[] 712 
00000 000 AY545738” AY545766" 
Qumarléb[] Qinghai [] J. 0. Liu[] 974 
0000 000 AY545736” | AY545767* 
Jomda[] Xizang 0 J. Q. Liu[] 1264 
00000 0000 000 O/CUS = AY545760* 
A. spinuliferd] Franch.[] R. Knuth Sangri[] Xizang [] J. Q. Liu[] 1138 
00000 0000 000 Q/CUS AY545747” — AYSASTIO" 
A. squarrosula Maxim. Xinghai[] Qinghai П J. Q. Liu[] 630 
000000 0000 000 e/cus - AY545772* 
A. tangulashanensis Y. C. Yang Zhidoi[] Qinghai [] J. Q. Liu[] 681 
& R. F. Huang 
00000 0000 000 Q/CUS AY545750 = 
A. tapete Maxim. Маф Qinghai [O J. Q. Liul] 627 
00000 0000 000 O/ROS _ Ау545749" AY545773 * 
A. wardi W. W. Smith Yushu[] Qinghai [] J. Q. Liu[] 884 
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000 ooo 00005 O00000%  Genbak0 00 
Taxon Locality Voucher! Distribution GenBank accession number! 
20] 
and habit 115 tmL-F 
00000 0000 000 Q/CUS AY545746 ü 
A. yargongensis Petitm. На Qinghai П J. Q. Liu[] 442 
0000 000 AY545743 * - 
Маф Qinghai 0 J. Q. Liu[] 523 
0000 000 AY545744” — AY545781 * 
Nanggén[] Qinghai [] J. Q. Liu[] 594 
0000 000 AYS45745' АҮ545782* 
Chindu[] Qinghai [O J. Q. Liu[]841 
00000 0000 000 Q/CUS _— AY545740' _ AY545768" 
A. zambalensid] Petitm.[] Hand.-Mazz. Maqén[] Qinghai [O J. Q. Liu[] 506 
0000 000 AY545741 ° АҮ545769* 
Zhidoi[] Qinghai [] J. Q. Liu[] 683 
00000 000 AY545742* = 
Qumarléb[] Qinghai [] J. Q. Liu[] 954 
A. alpind] L.[] Lam. O GenBank[] А /С05 AY275053 AY27496 
A. brevi] Hegetschw.[] Севан 0 GenBank[] А /CUS AY275049 AY274964 
A. chamaejasme Wulf. [] GenBank[] AJ ROS AF402555 AF402437 
ППП ПП A. delavayi Franch. [] GenBank[] @/ROS AY275101 AY274950 
А. heveticad] L.[] All. O СепВапк A /С05 AY275065 AY274981 
00000 A. hookeriana Klatt O GenBank[] O/ROS AY275096 AY274954 
A. lactea L. O GenBank{] А /С05 АҮ275069 AY274986 
UD UL A. lehmanniana Spreng. [] GenBank[] A/ROS AY275088 AY274948 
A. obtusifolia АП. 0 GenBank[] 4/ROS _ AY274972 
A. pubescens DC. 0 GenBank[] А /CUS AY275064 AY274978 
ПП ПП П] A. rigida Hand.-Mazz. [] GenBank[] O/ROS AY275095 AY274951 
А. sempervivoides Jacq. 0 GenBank] A /ROS AF402555 AF402436 
A. studiosorum A. Kress 0 GenBank[] A/ROS _ AY274946 
А. vandelli(] Тита Chiov. O GenBank{] A/CUS AY275060 AY274967 
A. villosa L. O GenBank{] А /С05 AY275084 AY274944 
A. wulfeniana Sieber ex Koch 0 GenBank[] А /CUS AY274962 AY275047 
00000 O СепВапк @/ROS AY275085 AY274952 
A. stenophylla{] Petitm.[] Hand.-Mazz. 
O00 Sect. Orthocaulon Hand.-Mazz. 
00000 0000 000 ФПО/АМАМ AY545757” AY545780 
A. erecta Maxim. Bangda[] Xizang 0 J. Q. Liu[] 1224 
UD Sect. Samuelia Schlechtd. 
00000 0 GenBank[] (О/РВЕ AY275105 AY274949 
А. axillari$] Franch.[] Franch. 
00000 0000 000 O/PRE — AY545748' АҮ545771° 
А. henryi Oliv. Jiulong[] Sichuan [] J. Q. Liu[] 670 
00000 O GenBank[] (О/РВЕ AY323705 - 
A. paxiana R. Knuth 
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000 000 00005 000002 — GenBankO OO © 
Taxon Locality Voucher! Distribution ^ GenBank accession number 
and habit ims ЊЕ 
00000 0000 000 ФПО/АМА AY545753” AY545776" 
Pomatosace filicula Maxim. Мад Qinghai П J. Q. Liu[] 628 
diri ан 0000 000 AY545752” AY545775" 
Chindu[] Qinghai 0 J. Q. Liu[] 649 
0000 000 AY545751” AY545774 * 
Zhidoi[] Qinghai 0 J. Q. Liu[] 686 
0000 000 AY545755” AY545777 ° 
Xinghai[] Qinghai [ J. Q. Liu[] 738 
0000 000 AY545754” _ AY545778 ` 
Madoil] Qinghai 0 J. Q. Liu[] 739 
Douglasia nivalis Lindl [] 2n = 380 0 GenBank[] ECUS AY275026 AY274991 
Vitaliana primuliflora Bertol [] 2n = 400 0 GenBank) А /CUS AY275050 AY274966 
0000 0 GenBank[] = AJ491428 = 
Bryocarpum himalaicum AJ491672 
Hook. f. & Thoms. 
Soldanella rugosa L. B. Zhang [] GenBank[] - AJ306360 = 
0000 П GenBank[] _ AJ491669 _ 
Stimpsonia chamaedryoides 
Wright ex A. Gray 
00000 0000 000 = AY545759 ° — AY545784 ° 
Primula gemmifera Batal. Nangqên[] Qinghai [] J. Q. Liu[] 613 
1ШПППППППППППППППППППППШ НР 2000 0 0 0000 00000 000 0 0 0000 
ОПОФООООООООООАПООФШОООООООАОО 000 0MOO00O000O0000 00 ПАМО0 0 
О 000000000000 О cusi 000 0 ОРВЕОО 0 0 0 козо 00000030 +00000000000 








10 Vouchers are preserved in НМР. 20 Distribution ОО Southeastern Qinghai-Tibetan Plateau and southwestern China] ФП Cen- 
tral Qinghai-Tibetan Plateau[] 20 Central Asia[] northern Asia and Caucasus[] АП Europe mountains] Ш North America. Нађи 
АММО annuals{] biennials or once seed set after growing for more than two years[] CUS[] cushion-like perennials[] PRE[] perennials[] 
ROS[] rosette perennials. 3[] * [] sequences are firstly reported in the present study. 


Tris-HCl [] O [I] pH 9.20 25 mmol/L KCI[] 1.5 mmol/L МС 5 pL dNTP 1 и] 8 mmol/L] 
00000 2.5 ull] 5 pmol/L Tag 0 1 pl) 3 Оту LOO 0 0 0 2 Ш 12.5 mmol/L 
00000 C4 тШ O 94 °С 40 «055 °С 20 s072 < 50 s00 0 32000 72 ©ПП 4 
min] 00000 PCR [|] O O O CASpure PCR Purification Kit] [] JUD U 
ПОООООШОООООООООООООО DYEnamic Dye Terminator Cycle Se- 
quencing Kit] Amersham Biosciences Corporation|] O 00000 РСКООШО О 04 < 15 
sLb0 °C 15 з060 °С 905] O 250 0 0 0 O Autoseg 96 Plated] Amersham Biosciences Corpo- 
ration |] O ШП MegaBACE 500 DNA Analysis ЗубетОДПО00000000000000О 
500 БР 0000000000 50%0 

TmL-F[] ITSI] 00000 СепвакО [] HL B] UO. D. D. D] UO D. D. 7. D. DO. U. UO LU. U 
000 GenBank T] O GenBank O O 0 0 U О DVitaliana [] O Douglasia DO D] B. UU. DIU. UI 
00000 GenBank [1 OOOO 10] 
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Fig. 1. Fifty percent Bayesian majority rule consensus tree of trnL-F sequence data. Numbers above internodes indicate the 
Bayesian posterior probabilities and below the MP bootstrap value based on 1000 replicates. Dashed lines denote that collapse in the 
strict consensus tree of parsimonious analysis. ANDR[] SAMU[] СНАМО AIZO[] and ORTH following each Androsace species are ab- 
breviations of sect. Androsace[] sect. Samuelia[] sect. Chamaejasme[] sect. Aizoidium[] and sect. Orthocaulon[] respectively. Habits 
of each species following section names are indicated as[] ANN[] annuals[] biennials or once seed set after growing for more than two 
yeas[] CUS[] cushion-like perennials[] РКЕП perennial] ROS[] rosette perennials. The distribution of each species is illustrated in the 
above left corner of figure. 











60 ПППППШПППППШПППШППППШППППППШПППППППП 489 





50%ПППППППППППППП so*#HO000000 00000000600 
ПШППШПППППППП 0000 "0 40000 2000 2000000000 
Bayesian ООООООООООООООООО Ооо О О 00 0000 00 
ШППШПШШППРПППППШПППППШППППППП1ПШПИПШПППП 
00000000 Ш0 0 ВауеһППППШПШПППППШПШППППППШППППП 


О Southeastern Qinghai-Tibetan Plateau and southeastern China Androsace zambalensis 506 CHAM CUS o 
100 A vis 594 
Q Central Qinghai-Tibetan Plateau 65 Еа CHAM cus e 
. 100 A. limprichtii CHAM ROS o 
AN Central Asia, northern Asia and Caucasus 91 A. zambalensis 683 CHAM CUS e 
100 7 " 
A European mountains 87 А. yargongensis 841 CHAM CUS o 
Bl North Ameri 100 91 А. chamaejasme CHAM ROS A 
9 menca 89 8! A. lehmanniana CHAM ROS @e 
100 A. studiosorum CHAM ROS A 
100 А. sempervivoides CHAM ROS A 
100 A. mariae 622 CHAM ROS eo 
100 94 A. wardii CHAM ROS Q 
90 
i И татав | СНАМ ROS г Ө) 
9 . mariae 
Т" a 9 CHAM ROS г Ө) 
7 50 А. татав 974 СНАМ ROS eo 
28 А. integra AIZO ANN/PRE® 
100 A. spinulifera CHAM ROS О 
100 А. stenophylla СНАМ ROS O 
99 Á. bisulca 1140 CHAM ROS O 
100 [ 65 A. bisulca 1125 CHAM ROS О 
100 Ç 
En i A. bisulca 1160 CHAM ROS O 
74 90 A. squarrosula CHAM CUS o 
100 63 A. alaschanica var. zadoensis CHAM CUS o 
100 A. tangulashanensis CHAM CUS e 
100 69 А. hookeriana CHAM ROS Q 
109 62 А. rigida CHAM ROS О 
Á. villosa CHAM CUS A 
А. delavayi CHAM ROS О 
А, erecta ORTH ANN ФО 
А. alpina CHAM cus A 
Ст А. hevetica СНАМ CUS A 
92... A. obtusifolia CHAM ROS A 
56 A. lactea CHAM cus A 
Douglasia nivalis CUS m 
Ted NE 98... [ев CHAM с А 
68 59 A. wulfeniana CHAM CUS A 
100 A. pubescens CHAM CUS A 
100 190 Vitaliana primuliflora CUS A 
100 A. vandellii CHAM CUS A 
A, septentrionalis ANDR ANN AA 
оо T EE Pomatosace. filicula 686 ANN eo 
P. filicula 739 AN ФО 
100 P. filicula 628 AN ФО 
5 100 Р. filicula 738 AN ФО 
86 P. filicula 649 ANN eo 
A. henryi SAMU PRE O 
A. axillaris SAMU PRE О 
А. maxima ANDR ANN A 
4. filiformis ANDR AN A 





Primula gemmifera 





490 uU 00000: 42 





ügggaaggaggauaggugma 80 000020 00000000 utu 
0 O UDouglasia nivalis П Vitaliana primuliflora [] LI] 00000-00000 00000 
ОФ оОоООоОоОоОООооооошоооооооооооооо О0О О LDougla- 
sia nivalis П Vitaliana primuliflora DU OO Ш 00000000 20000000000 
ПШШШПШПШПППППШПШПШПППППШППППШПППППППШППППШПШПШПП 
00 

msputuutiull 554 БРПОООООООШОО 0 00 3200 Ш 00 
000000 5 00000 12090 000 0 0ОС/ = 0. 5580! = 0.81300 0 = 12240 
Bayesian 50%0 [] HH ПО ООО 22D udi ioo Dag d Bg Bg B OU UO UU D U 
HOUUMOUUUUU S%ILULUDUDUDUD28UUUUUUDIUDUUUUUUUUunnmu 
umnpmuumumg ПППП”ПППППППППППШШПППП ПП Bayesian 0 0 
0000 USoldanella тироза Stimpsonia chamaedryoides JO D] D] 0. d] D] B. D] L7. 0000" 
OOOOOUU 0000 "О 0150 0000000 ттроооооооооооо 
HOUOUUU 3ЗППППППШПШПШПШППППШППППШШПШПШППППШППППППП 
OOOO 08000 "0000000 serragumuuagguuar 00000-00 
EL ELSE ЕЦ 

ШШППШПППШПШШППШШ D Ws+mLFOU0U00U000 45000000 
1645 БР 00000000000000 4100000000000000000 
ПООООООООООООООООО ПО 0 РРОЦООООООООООШО 
ООООООООШ Р=0.80 >0.05Ш ППППППППП 400000 400 
000 0С/ =0.6620 RI = 0.87000 П = 15120 Bayesian 50%0 0000000 30000 
ШПШПШПШПППШПППШПППШППППШППППППППШШППШППППППППППППП 
О0ОоОоОо" ОООО "ВОООООЙОООО "0" ППППП-ПППП”ППППП 
ШППШПШШПШШПППШПППШПППШППППШППППППППППППШПШПШПППШПШПП 
OOOOOM ПППП”ПППШППППППППППППППППП 
2.2 000000000 

O O Modeltest 3.04] Posada & Crandall] 1998[T] 0 0 0 500000000 СТК + 
ctn 000000 ШОООООООО О among-site rate variation] gamma O [T]. gamma 
distribution [T] 0 0 O alpha= 1.033500 0 0 0 ШП substitution rate[] &] all = 1.21190 R] bO 
= 2.2055) c= 1.587708 40 = 0.450008) e[] = 3.668400 O 00 0 ПП base frequency А 
= 0.2142[С = 0.272406 = 0.26950T = 0.24390 00 0000000000000000 
0000000000000000 Ш - 2А = 2] 5387. 31218 – 5330.569280 = 
113.48580 ај = 510 P <0.005 A0 0000000000000000 0 0 0( 0 00 0 0 
iDpBugugugmmpuaumugpguulbuituuuututututututtu iu t 
00000000 7.83x10^?s s^* a T] substitution per site per уе | Sang et al. [] 199500 
OM ПППП”ППППППППППППП ПШ 18.02 Му] million years ago] ПП 
ООШОШО О мка 0 000000000000 40 0 00 0 0 0 00 0 00 00 
П0000000 овоо оовоооовооооооо 2000000 





+ 
© 
= 


60 ООООоООоОоОоОоОоОоОооошоооооооооооооо 





CO Southeastern Qinghai-Tibetan Plateau and southeastern China 








100 Androsace zambalensis 954 CHAM CUS 

@ Central Qinghai-Tibetan Plateau A. yargongensis 442 CHAM CUS 
A. yargongensis 523 CHAM cus 

AN Central Asia, northern Asia and Caucasus A. yargongensis 841 CHAM  CUS 
A European mountains 4. yargongensis 594 CHAM CUS 
А. zambalensis 506 CHAM CUS 

m North America zambalensis 683 CHAM CUS 
chamaejasme CHAM ROS 

lehmanniana CHAM ROS 

bisulca 1140 CHAM ROS 

bisulca 1160 CHAM ROS 

bisulca 1125 CHAM ROS 

sempervivoides CHAM ROS 

squarrosula CHAM CUS 

alaschanica var. zadoensis CHAM ROS 

limprichtii CHAM ROS 

mariae 1264 CHAM ROS 

wardii CHAM ROS 

mariae 622 CHAM ROS 


mariae 974 CHAM ROS 
mariae 712 CHAM ROS 
integra AIZO ANN/P 
tapete CHAM CUS 
delavayi CHAM ROS 
stenophylla CHAM ROS 
villosa CHAM CUS 
hookeriana CHAM ROS 
rigida CHAM ROS 
erecta ORTH ANN 
maxima ANDR ANN 
filiformis ANDR ANN 
hevetica CHAM cus 
A. pubescens CHAM CUS 
А. lactea CHAM ROS 
Douglasia nivalis CUS 
A. vandellii CHAM CUS 
A. alpina CHAM CUS 
A. wulfeniana CHAM cus 
A. brevis CHAM CUS 
Vitaliana primuliflora CUS 
A. septentrionalis ANDR ANN 
Pomatosace filicula 739 ANN 
P. filicula 649 ANN 
P. filicula 686 ANN 
P. filicula 738 ANN 
P. filicula 628 ANN 
A. paxiana SAMU PRE 
A. axillaris SAMU PRE 
A. henryi SAMU PRE 


Bryocarpum himalaicum 


~RRBARRRARARARARRARAR вв RD 


00099990 <<. »»И»»»>>аоовеееваваовеве> ооо>Ревевеве 


Soldanella rugosa 
Зитрхота chamaedryoides 


Primula gemmifera 


02 WSOO0 Bayesin 50%0 00000 000000 0 


Fig. 2. Fifty percent Bayesian majority rule consensus tree of ITS sequence data. The remaining explanations are the same as in 


Fig. 1. 
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Fig. 3. Fifty percent Bayesian majority rule consensus tree of trnL-F + ITS sequence data. The remaining explanations are the same 
as in Fig. 1. 
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Fig. 4. The crudely estimated divergence time of main branches according to the NPRS] nonparametric rate smoothing[] tree cali- 
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